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Resonant magnetic excitations are widely recognized as hallmarks of unconventional superconductiv- 
ity in copper oxides/'^ iron pnictides,^ and heavy-fermion compounds.^"^ Numerous model calcula- 
tions have related these modes to the microscopic properties of the pair wave function/"^^ but the 
mechanisms underlying their formation are still debated.^ Here we report the discovery of a simi- 
lar resonant mode in the non-superconducting, antiferromagnetically ordered heavy-fermion metal 
CeB5.^^"^^ Unlike conventional magnons, the mode is non-dispersive, and its intensity is sharply con- 
centrated around a wave vector separate from those characterizing the antiferromagnetic order. The 
magnetic intensity distribution rather suggests that the mode is associated with a coexisting order pa- 
rameter of the unusual antiferro-quadrupolar phase of CeB^,^^^^ which has long remained "hidden" 
to the neutron-scattering probes. ^^'^^ The mode energy increases continuously below the onset tem- 
perature for antiferromagnetism, T^, in parallel to the opening of a nearly isotropic spin gap through- 
out the Brillouin zone. These attributes bear strong similarity to those of the resonant modes observed 
in unconventional superconductors^"^ below their critical temperatures, T^.. This unexpected common- 
ality between the two disparate ground states indicates the dominance of itinerant spin dynamics in 
the ordered low-temperature phases of CeB5 and throws new light on the interplay between anti- 
ferromagnetism, superconductivity, and "hidden" order parameters in correlated-electron materials. 



In copper-oxide, iron-pnictide, and heavy-fermion supercon- 
ductors, antiferromagnetic (AFM) order can be induced by 
changes of pressure, doping level, or magnetic field. The- 
ories of magnetically mediated Cooper pairing inspired by 
this observation have received strong support from neutron 
scattering experiments that revealed "resonant" magnetic 
excitations in the superconducting states of all of these 
compounds. The resonant modes share a common set 
of characteristics. They are sharp in both energy and mo- 
mentum, unlike the incoherent particle-hole excitations in a 
metal, and are located below the superconducting gap, 2 A. 
Since model calculations indicate that the spectral weight 
and momentum-space structure of the resonant modes con- 
tain essential information about the pairing interaction, ^'^'^ 
pairing symmetry, and condensation energy^^ of the su- 
perconducting state, research on their origin remains at 
the frontier of the quest for a comprehensive description of 
unconventional superconductivity. 

Here we report inelastic neutron scattering (INS) exper- 
iments that reveal a resonant mode with identical charac- 
teristics in the non-superconducting heavy-fermion metal 
CeBg, which has been intensely investigated because of its in- 
triguing phase behavior^ on the background of an excep- 
tionally simple crystal structure and chemical composition. 
The electronic states near the Fermi level of CeBg are com- 
posed of localized cerium-4/ ^ levels hybridized with itin- 
erant cerium-Sd and boron-2p electrons, ^^'^^ closely analo- 
gous to those of the superconductors CeCoIUs and CeCu2Si2, 
where prominent resonant modes have been reported. "^'^ 
In contrast to these superconducting compounds, however. 



CeBg exhibits AFM order below = 2.3 K.^^ The AFM tran- 
sition is preceded by another phase transition at Tq = 3.2 K, 
whose order parameter has long remained "hidden" to stan- 
dard experimental probes such as neutron diffraction.^^'^^ 
Theoretical work^^'^^ has attributed both transitions to in- 
teractions between the multipolar moments of the Ce-4/ 
electrons mediated by the itinerant conduction electrons. 
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Figure 1 | Bragg peaks associated with the AFQ and AFM or- 
der parameters, a, Comparison of the elastic signals measured at 
Qafq = ^( III) and Q^fm = |0) wave vectors at T = 1.6K. b, 
Temperature dependence of the peak intensities that vanish at Tq 
and T^, respectively. The solid lines are guides to the eyes. 
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Figure 2 | Unprocessed INS data, a, Inelastic signal along the (HHL) plane in reciprocal space, measured in the low-temperature AFM 
phase (T = 1.6 K) at 0.5 meV energy transfer. The labels represent different high-symmetry points, and the dashed lines mark Brillouin 
zone boundaries, b, c. The same at lower (0.25 meV) and higher (0.9 meV) energies, respectively, d. Energy dependence of the INS signal 
in the AFQ (empty symbols) and AFM (filled symbols) states, measured at the R(^^^) point, where the intensity is maximized. Exciton 
peaks are marked by arrows. The horizontal black bar inside the peak shows the experimental energy resolution, the hatched region 
represents background intensity, and the lines are guides to the eyes. The inset shows constant-energy scans across the peak maximum 
(Hco = 0.5 meV) along the MRM^ line at different temperatures; the dashed line is the background level suggested by the intensity minima 
at M(||0) and M'(||l) within the spin gap (Hco = 0.25 meV). e. Low-temperature (T = 1.6K) energy scans for theX, M, T, and A 
points. The location of the scattering plane and the naming of high-symmetry points in the cubic Brillouin zone are shown in the inset. 



which break the large ground-state degeneracy of the cerium 
ions in their cubic crystal field and stabilize an antiferro- 
quadrupolar (AFQ) order. 

Recent transport, thermodynamic, and neutron diffrac- 
tion experiments^^'^"^"^^ have revealed deviations from pre- 
dictions of this "interacting multipole" model of magnetism 
in CeBg. In particular, w^eak Bragg intensity has been de- 
tected by spin-polarized neutron scattering at the AFQ wave 
vector, although the concept of quadrupolar ordering for- 
bids such an intensity in the absence of external magnetic 
fields,^^'^^ so that it is "hidden" to the elastic neutron scat- 
tering. This observation has been ascribed to the magnetic 
response of the itinerant Ce-Sd electrons. Our discovery 
of a magnetic resonant mode for T < confirms a more 
pronounced influence of itinerant electrons on the magnetic 
properties of CeBg than anticipated by the prevailing theory, 
and establishes a surprising correspondence to unconven- 
tional superconductivity. 

We start the presentation of our data by comparing the 
elastic neutron scattering intensities measured at the pre- 
viously identified AFM and AFQ Bragg peaks at Qafm — 
^0) and at Q^fq = respectively, in Fig. 1. From 

the different onset temperatures of the elastic signals, it is 
clear that they represent distinct order parameters. The 
elastic intensity at Qafq is more than an order of magni- 
tude weaker than the AFM Bragg intensity, and exhibits only 
a w^eak suppression below^ T^, indicating that both order 
parameters coexist in the AFM phase. 

To give an overview of the inelastic-scattering signal, we 
present in Fig. 2a-c unprocessed INS intensity maps of the 
(HHL) plane in reciprocal space, measured in the AFM low- 
temperature phase (r = 1.6K) at fixed energy transfers of 
Hco = 0.25, 0.5, and 0.9 meV An intense peak centered at 
Qafq dominates the spectrum at Hco = 0.5 meV It can be 



observed at the two equivalent R points (HI) and -|) 
covered by the measurement, w^hich are additionally con- 
nected by weaker streaks of intensity that surround the M 
point. Intensity maps at both lower (Fig. 2b) and higher 
(Fig. 2c) energies reveal greatly reduced intensities. The 
strong concentration of magnetic spectral weight around 
isolated points in energy-momentum space is markedly dif- 
ferent from the continuously dispersive spin waves in con- 
ventional antiferromagnets, but closely similar to the be- 
havior of the resonant modes observed in copper-oxide, 
iron-pnictide, and heavy-fermion superconductors below 
their respective transition temperatures.^"^ Moreover, the 
mode is not centered at the S or 17 points characterizing the 
AFM order, where the magnetic intensity is much weaker, 
but at the R point that reflects the "hidden" AFQ phase. 

At T = 1.6 K, the resonant mode is sharp in energy 
(Fig. 2d) and is characterized by the intrinsic (resolution- 
corrected) full width at half maximum (FWHM) of 
~ 0.2 meV It is followed by a weaker peak centered around 
1.1 mey with the FWHM of ~0.6meV For T > T^, the 
0.5 meV peak transforms into a quasielastic response^"^'^^ 
that is typical for the paramagnetic state of heavy-fermion 
metals. We can therefore associate the formation of 
the resonant peak with the opening of a spin gap in the 
low-energy part of the quasielastic spectrum upon cooling 
and the consequent spectral-weight transfer from below 
~ 0.35 meV to higher energies. In contrast, the broader 
peak at 1 . 1 meV is nearly insensitive to and is possibly 
due to a low-energy crystal-field excitation similar to those 
previously studied in the AFQ phase. ^^'^^ 

The detailed temperature dependence of the resonant 
mode is illustrated by Fig. 3, where we plot the imaginary 
part of the dynamical spin susceptibility, ;^^^(Qafq^ ob- 
tained by correcting the background-subtracted INS signal 




Figure 3 | Temperature dependence, a, Imaginary part of the dynamic spin susceptibility, j^^CQafq? various temperatures, 

illustrating the gradual development of the spin-exciton peak below the AFM transition. The lines are visual guides, b, Temperature 
dependence of X^XQafq^ ^) the peak maximum (Hco = 0.5 meV) and within the spin gap (Hco = 0.25 meV). The inset shows the 
exciton energy versus temperature, determined from the position of the peak maximum, Hco^^^, in panel a. c. Comparison of the 
normalized resonance energies, hco^^^/k^T^ (or TJ, plotted vs reduced temperature, T/T^ (or T^), in CeBg, in several unconventional 
superconductors,^'^'^ and in the low-temperature phase of CeRu2Alio-^° The lines are guides to the eyes. 



for the thermal population factor. In Fig. 3a, one sees that 
the mode gradually broadens and shifts towards lower en- 
ergies upon warming, following the order-parameter-like 
trend shown in the inset of Fig. 3b. This behavior is again 
strikingly similar to the one of the resonant mode in un- 
conventional superconductors (Fig. 3c). For T > T^, the 
quasielastic signal continues to diminish in intensity with 
increasing temperature (Fig. 3b) . Constant-energy scans 
in the inset of Fig. 2d demonstrate that it always remains 
peaked at the commensurate wave vector Qafq. where the 
resonant mode develops below T^. 

Constant-Q scans at various high-symmetry points 
(Fig. 2e) further show a weakly energy dependent mag- 
netic response throughout the Brillouin zone. For T < T^, a 
spin gap of magnitude ~ 0.5 meV opens up nearly isotrop- 
ically in momentum space, including the M point where 
the magnetic INS signal is minimal, and the R point where 
the resonant mode develops in the AFM state. The appear- 
ance of the spin gap is consistent with the depletion in the 
electronic density of states below T^, previously observed 
by point-contact spectroscopy."^^ The comparison in Fig. 3b 
of the temperature scans measured at the peak maximum 
(0.5 meV) and below it (0.25 meV) emphasizes that the 
spectral-weight gain at the exciton energy is compensated 
by its depletion within the spin-gap region. Both effects 
show a simultaneous sharp onset at T^, closely similar to 
the behavior of the resonant mode and the spin gap in an un- 
conventional superconductor.^ An analogous T -evolution is 
observed at the X(00|) point (see Fig. S2 in the Supplemen- 
tary Information), where a quasielastic signal of comparable 
intensity is found in the AFQ state. The redistribution of 
spectral weight at this point leads to an intensity enhance- 
ment below at hco ^0.9 meV, but the resulting feature is 
much broader than the one at the R point. 

Although the resonant mode is most intense at R( ), 
Fig. 2a suggests the presence of weaker tails extending well 



into the Brillouin zone. This observation is substantiated 
by Fig. 4, where we show momentum scans along the ATA^ 
direction that cut the two lobes of intensity extending from 
R towards the M point (see also Fig. S2b in the Supplemen- 
tary Information) . The comparison of the low-temperature 
signals measured at 0.5 and 0.25 meV together with the 
energy scan at the A point (Fig. 2e) let us conclude that 
the mode exhibits at most a weak dispersion, so that even 
away from the R point its intensity remains maximal near 
^^max ^0-5 meV 

While the quasielastic intensity in the paramagnetic state 
of CeBg is similar to the one observed in other heavy- 
fermion compounds,^'^ the behavior of the resonant mode 
we have discovered for T < is quite different from that of 
low-energy magnons in other heavy-fermion antiferromag- 
nets including the isostructural PrB5"^^ as well as Celug"^^ 
and CePd2Si2."^'^ These excitations exhibit the conventional 
magnon dispersion with minimum energy and maximum 
intensity at the AFM ordering wave vector. The resonant 
mode in CeB5, on the other hand, is much more intense, at 
most weakly dispersive, and centered at Qafq rather than 
Qafm- The behavior we have observed disagrees with calcu- 
lations in the framework of the interacting-multipole model 
of CeBg, which predict magnon-like low-energy modes,^^ in 
addition to more weakly dispersive higher-energy crystal- 
field excitations.^^ 

We have already pointed out the striking similarity of the 
resonant mode in CeBg to the low-energy spin excitations 
in unconventional superconductors. While this analogy has 
not been anticipated as far as we are aware, we can under- 
stand it qualitatively in the framework of a two-component 
model that includes strong interactions of the localized mag- 
netic moments with the itinerant electrons.^ The formation 
of the AFM state accounts for the opening of a gap in the 
electronic density of states near the Fermi energy, which is 
similar to that found in spin-density-wave (SDW) systems. 



3 




0.2 0.4 0.6 0.8 

Hin (HH 1/4) 



Figure 4 | Anisotropic shape of the resonant mode in momen- 
tum space. Low-temperature Q-scans along the ATA! direction 
cross the two lobes of the resonant intensity, resulting in two in- 
commensurate peaks. The intensity is maximized at the resonance 
energy (fico = 0.5 meV) and is much reduced within the spin gap 
(fico = 0.25meV). 

As a consequence, it leads to the nearly isotropic spin gap in 
the continuous magnetic spectrum originating from strongly 
hybridized crystal-field excitations of the Ce 4/ -electrons 
and the itinerant conduction electrons. The resonant mode 
can then be regarded as a collective mode below the onset 
of the particle-hole continuum, in close analogy to current 
theoretical models of the resonant modes in the supercon- 
ductors. We note that this scenario is quite different from 
the standard interacting-multipole model of CeB5, accord- 
ing to which the magnetic response in the AFM phase is 
dominated by the magnetic dipole moments of the localized 
Ce electrons, but qualitatively consistent with other recent 
observations that emphasized the importance of itinerant 
magnetism in this system. 

While there is general agreement on the collective-mode 
nature of the resonant modes in unconventional supercon- 
ductors, the mechanisms leading to their formation are still 
controversial. ^'^^ It is interesting to consider our data on 
CeBg in the light of two classes of theories that describe 
the mode as a magnon-like excitation of localized electrons 
close to a magnetically ordered state, and as an excitonic 
bound state of itinerant electrons, respectively. A recent 
calculation in the framework of the former approach^ has 
reproduced the temperature evolution of the magnetic re- 
sponse of superconducting CeCoIns from an overdamped 
mode above to a resonant collective mode below by 
assuming that the superconducting gap removes damping 
channels for a low-energy Ce crystal-field excitation. An 
analogous model in which the superconducting gap is re- 
placed by a SDW-type gap may explain the resonant mode 
in CeBg. On the one hand, the appearance of the reso- 
nant mode at the R point is consistent with the presence 
of AFQ order, which is expected to support a soft multi- 
polar mode at this wave vector. On the other hand, the 



close proximity of the resonant mode to the gap energy, 
which would have to be regarded as a coincidence in this ap- 
proach, is naturally explained in the framework of a scenario 
in which the mode is described as a weakly bound triplet 
exciton. The complex anisotropic shape of the collective 
mode in CeBg also suggests an influence of the fermiology 
of the conduction-electron system, in analogy to itinerant- 
electron descriptions of the resonant mode in CeCoIns^^ and 
the "hourglass" shape of the resonant mode in the copper 
oxides.^ While such an influence is generally expected for 
an exciton within an itinerant-electron description, it will 
be interesting to explore theoretically the impact of the co- 
herence factors on the expression for x"{S)., co) in the AFM 
ordered state. This is important because the relationship 
between the superconducting coherence factors and the res- 
onant mode has figured prominently in the ongoing debate 
about the Cooper pairing symmetry in iron pnictides"^'^'^ and 
heavy-fermion compounds^'^^. Answers to these questions 
will have to await calculations that take full account of the 
strong coupling between the localized Ce crystal-field states 
and the itinerant electrons, qualitatively beyond the current 
interacting-multipole model of CeBg. 

We end our discussion by noting that the resonant mode 
in CeB5 apparently results from the coexistence of the AFM 
state with the unusual AFQ order that had remained "hid- 
den" to experimental probes in zero magnetic field until the 
recent advent of resonant x-ray diffraction methods. 
This unique situation explains its unusual behavior in com- 
parison to other heavy-fermion antiferromagnets. Presum- 
ably as a consequence of AFQ -order fluctuations, the dynam- 
ical susceptibility at the R-point of CeBg remains substantial 
even for T > T^, where the resonant mode is overdamped, 
and only decreases gradually upon further heating. This be- 
havior is reminiscent of the gradual onset of magnetic inten- 
sity in the "pseudogap" state of the underdoped cuprates,^ 
for which various "hidden" order parameters have also been 
proposed. The mode we have observed may also be re- 
lated to the higher-energy excitonic modes that have been 
reported in other materials with thus-far unidentified "hid- 
den" order parameters, such as CeRu2Alio.^^ and in small- 
gap "Kondo insulators" SmBg^^ and YbBi2,^^ whose ground 
states remain poorly understood. Since the "hidden" order 
in CeB5 has actually been extensively characterized, the 
comprehensive description of the magnetic dynamics of this 
compound we have initiated here will have model character 
for a broad class of correlated-electron materials. 

Methods. 

Two single-crystalline rods of CeB5 with a total mass of 8 g were grown by 
the floating-zone method from a 99.6 at. % isotope-enriched ^^B powder, 
in order to minimize the absorption of neutrons by the ^^B isotope. 

The data were collected using the cold triple-axis IN14 spectrometer 
(ILL, Grenoble, France) with a pyrolytic graphite (PG) monochromator. 
The sample was mounted into a standard cryostat with the (110) and (001) 
directions in the scattering plane. To acquire the overview maps of the 
reciprocal space (Fig. 2a-c), we used the Flatcone detector, equipped with 
Si (111) analyzers that select scattered neutrons with the finite momentum 
kf = 1.40 A~^. In this configuration, a cold beryllium filter was installed 
in the path of the incident beam to eliminate the contamination from 
higher-order neutrons. As a consequence, the accessible energy- transfer 
range was restricted to 0.9 meV by the Be-filter cutoff energy. 
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The rest of the measurements were done in the conventional triple-axis 
configuration with a PG analyzer and the Be filter installed on k^. The 
same fixed value of fcf = 1.40 was selected. The data were corrected 
for the energy-dependent fraction of higher-order neutrons on the monitor. 
The imaginary part of the dynamical spin susceptibility, j^^(Q, co), was 
obtained from the background-subtracted INS intensity by the fluctuation- 
dissipation relation x'XQ, = (1 - e~^'^^^^'^)SiQ, co), where S(Q, co) is 
the scattering function. The error bars in all figures correspond to one 
standard deviation of the count rate. 

Throughout the paper, we index the reciprocal-lattice vectors on the 
simple-cubic unit cell (space group Pm3m, lattice constant a = 4. 1367 A) 
and follow the conventional notation in labeling the high-symmetry points 
(Fig. 2e, inset) . The wave-vector coordinates are given in reciprocal lattice 
units (1 r.l.u. = 2% /a). 
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Figure SI | Temperature dependence of the inelastic intensity at Kco = 0.5 meV The exciton peak at the point is significantly 

suppressed at T = 2.6 K > and nearly vanishes upon further warming into the normal state (T = 10 K), in agreement with Fig. 3. All 
panels share the same color scale. 



Here we present additional INS data to substantiate the temperature dependence of the new exciton mode and its shape 
in the reciprocal space. Additional maps presented in Fig. SI above demonstrate that the strong intensity at the R point 
(panel a) is sharply suppressed above (panel b), where it transforms into the quasielastic response that continues to 
diminish with further temperature increase (panel c) . 

In Fig. S2a, we directly compare energy scans (raw INS intensity) at the X and A points, measured in the low-temperature 
AFM state (T = 1.6K) and in the AFQ phase above (T = 2.6K). The low- temperature datasets are identical to those 
shown in Fig. 2e. At both points, the transfer of the low-energy spectral weight into a pile-up peak above the spin gap can 
be observed. At the A point, this peak is centered at ~ 0.5 meV and represents the "tail" of the resonant excitation centered 
at R, which demonstrates the absence of any notable dispersion of the exciton mode. 

In Fig. S2b, we show additional constant-energy scans along the (HHO) direction, centered at the M point. Despite 
the fact that this scan is shifted from the R point by half the Brillouin zone size, the peak intensity is still maximized at 
0.5 meV, with a deep minimum at the M point. This indicates that the two nearly non-dispersing lobes of intensity that 
extend from the exciton peak (see Fig. 2a and Fig. 4) close into an elongated donut-shaped structure in the (H, L,co) space, 
which surrounds the M point and is restricted to a very narrow energy range. At a lower energy of 0.25 mey the intensity 
is partially suppressed due to the presence of the spin gap, whereas at higher energies the Q-dependence of the intensity 
rapidly vanishes, as illustrated here by the 0.8 meV data. 
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Figure S2: a, Comparison of the 
energy scans measured at the X 
and A points above and below T^. 
b, Comparison of low-temperature 
momentum scans along the (HHO) 
direction, centered at the M point, 
for various energies. The scan at 
0.5 meV crosses the two intensity 
lobes that connect the resonant 
peaks at equivalent R points. 
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